Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN M?eX style file vl.4) 



SCUBA observations of galaxies with metallicity 
measurements: a new method for determining the relation 
between submm luminosity and dust mass 



A. James, L. Dunne, S. Eales and M. G. Edmunds 

Department of Physics & Astronomy, University of Wales, P.O. Box 913, Cardiff CF24 3YB, UK 



(N 
O 
O 

< 

O 

m 



> 
in 

o 

(N 
O 

6 



1 February 2008 



ABSTRACT 

Using a new technique we have determined a value for the constant of proportionahty 
between submilhmetre (submm) emission and dust mass, the dust mass absorption 
coefficient {nd) at 850/im. Our method has an advantage over previous methods in 
that we avoid assumptions about the properties of dust in the interstehar medium. 
Our only assumption is that the fraction of metals incorporated in the dust (e) in 
galaxies is a universal constant. To implement our method we require objects that 
have submillimetre and far-infrared (FIR) flux measurements as well as gas mass and 
metallicity estimates. We present data for all the galaxies with suitable measurements, 
including new submm maps for five galaxies. We find Kgso = 0.07 ± Q.Q2m?'kg^^ . 

We have also been able to use our sample to investigate our assumption that £ is a 
universal constant. We find no evidence that e is different for dwarf and giant galaxies 
and show that the scatter in e from galaxy to galaxy is apparently quite small. 

Key words: dust,extinction-galaxies:dwarf-galaxies:abundances-galaxies:evolution 



1 INTRODUCTION 



Measurements of the depletion of heavy elements, relative 
to the solar abundances, along the line of sight to the stars 
show that about 50% of the metals in the Galaxy is bound 
up in dust grains (Whittet 1992). This raises the questions: 
Is this fraction the same for every galaxy? What are the 
processes that so efficiently lock metals into dust grains? 

Although the presence of dust manifests itself in a mul- 
titude of ways - extinction, reddening, the polarisation of 
starlight, to name but three - the best way, in principle, to 
estimate the mass of dust in a galaxy is from the submm 
flux density. The advantages of this are that (a) the emis- 
sion from dust is optically thin and (b) the emission depends 
mainly on the mass of dust and is only weakly dependent 
on the temperature of the dust. The practical obstacle to 
making the full use of this technique has been the diffi- 
culty of estimating the constant of proportionality between 
the emission and the mass of dust, the mass-extinction or 
the mass-absorption coefficient. There are two methods that 
have been used to estimate this constant. 

The most commonly used technique to calibrate this 
constant is by observations of Galactic reflection nebulae 
(Hildebrand 1983; Casey et al. 1991). By making optical 
and far-infrared/submillimetre observations, one can esti- 
mate Av/r submm lu a reasonably straightforward way. If 
one assumes a value of Nh /Av and a gas-to-dust ratio, one 



can estimate the mass-absorption coefficient, k^. For obser- 
vations at long wavelengths such as with the SCUBA submil- 
limetre camera (Holland et al. 1999), which operates mostly 
at 850/xm, there is also the practical problem that this tech- 
nique has only been applied at relatively short (A < 400/im) 
submillimetre wavelengths, which means that the coefficient 
has to be extrapolated to longer wavelengths. Since the dust 
emissivity index, /?, (kj oc !/'') is uncertain by a factor of two, 
there is an immediate uncertainty of a factor of two in the 
long-wavelength value of Kd on top of all the other uncer- 
tainties. This problem, of course, could be solved by SCUBA 
observations of reflection nebulae, but there are additional 
problems with this technique. First, there is evidence that 
the properties of dust vary with environment (e.g. Cardelli 
et al. 1996), and so the dust in a reflection nebula may not be 
representative of a galaxy as a whole. Second, the technique 
relies on the assumption of spherical symmetry for the re- 
flection nebulae and on the values for properties that are not 
themselves simple to measure, such as the dust-to-gas ratio 
and Nh /Av (although there are variations in this technique 
that do not require these, there are other assumptions that 
have to be made, such as the typical sizes and densities of 
dust grains). Hughes et al. (1993) estimate that the factor 
of uncertainty in Kd{v) / [Mg / Md) at SOO/xm is ^ 25. 

Alton (2000) employed Hildebrand's (1983) technique 
using SCUBA observations of an external galaxy to derive 
Av/ffi^Q, essentially using a galaxy rather than a reflection 
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nebula to calibrate Kd- This method obviously avoids the 
problems of extrapolation to longer wavelengths and the fact 
that nebular dust may not be the same as dust in the rest 
of the interstellar medium but still requires assumptions to 
be made about the physical properties of the dust. 

The other technique that has been used to estimate Kd 
is more loosely linked to submillimetre observations. In this 
technique one constructs a model for dust (its chemical com- 
position, distribution of grain sizes etc.) that agrees with all 
the available data, for example the shape of the optical ex- 
tinction curve, and then uses the model to estimate the value 
of the submillimetre mass-absorption coefficient (Draine & 
Lee 1984; Hughes et al. 1993). 

In this paper we suggest an alternative way of estimat- 
ing the mass-absorption coefficient based on the global prop- 
erties of galaxies rather than on individual properties within 
the Galaxy. This technique is based on the assumption that 
the fraction of metals within the interstellar medium of a 
galaxy that is bound up in dust is a constant. There is 
some evidence for this from the correlation between dust-to- 
gas ratio and metallicity in nearby galaxies (Issa, MacLaren 
& Wolfendale 1990). The high depletion of many elements 
within our own ISM (Whittet 1992) also suggests that the 
mechanism forming dust is efficient - that is, if there are met- 
als present, they are efficiently incorporated into dust grains. 
Finally, there is evidence in the Galaxy that although the 
dust extinction curve is quite variable from place to place, 
the gas phase carbon and oxygen abundances are remark- 
ably constant over a wide range of ISM density (Cardelli et 
al. 1996; Meyer et al. 1998), implying (as these elements are 
the most important constituents of dust) that the fraction 
of metals in dust is a constant. In contrast to this evidence, 
Lisenfcld & Fcrrara (1998) found a non-linear relationship 
between metallicity and dust-to-gas ratio for 28 dwarf irreg- 
ular galaxies, which suggests that e, the fraction of metals 
locked up in dust, is not the same for all types of galaxy. In 
this paper, we will present new data that shows that Liscn- 
feld & Ferrara's technique missed cool dust, and when this is 
taken into account e does appear to be remarkably constant 
from galaxy to galaxy. 

To apply the technique we need a galaxy which has sub- 
millimetre and far-infrared ffux measurements, a measure of 
its total gass mass, and a measure of its metallicity. 

The dust mass of the galaxy is then given by two equa- 
tions: 

Md = MgX Z xex f (1) 

where Mg is the mass of gas in the galaxy, Z is the metallicity 

relative to solar, e is the ratio of the mass of metals in the 
dust to the total mass of metals and / is the ratio of the mass 
of metals to the mass of gas for gas with solar metallicity, 
we estimate / = 0.019, and; 

Md= (2) 

K860 X Bs50{l^,T) 

in which 5*850 is the ffux density at 850/xm, D is the dis- 
tance of the galaxy, Bs5o{i^,T) is the value of the Planck 
function at 850/^m, T is the dust temperature and resso is 
the mass-absorption coefficient at 850/im. We then set these 
two equations equal to one another, and as long as we know 
e, we can estimate Kd- 

The arrangement of this paper is as follows. In section 2 



we present the results of a search through the archive of the 
James Clerk Maxwell Telescope for submm data for galax- 
ies with known metal abundances. In section 3 we derive 
a value for e based on observations of the local interstellar 
medium. In section 4 we apply the technique to estimate Kd 
and investigate whether e is a universal constant. We assume 
throughout that Ho = 75kms~^ Mpc~^ . 



2 DATA REDUCTION & RESULTS 

We looked through the JCMT* archive for all the objects for 
which there are both abundance measurements and obser- 
vations with the SCUBA camera (Holland 1999). Table 1 is 
a list of galaxies in the archive for which there are published 
metallicity measurements but no published submm results. 

The data for these ffve objects was reduced in the usual 
way using the SURF (Jenness & Lightfoot 1998) software 
package. Images were first flat-fielded to remove inhomo- 
geneity in bolometer sensitivity. Corrections were then made 
for atmospheric absorption using the opacities (rsso) derived 
from skydip measurements. Noisy bolometers were masked 
and large spikes in the bolometer time-stream removed by 
applying a 4-a clip. Residual sky noise was removed with 
the SURF task REMSKY. FoUowing removal of the sky 
noise, the bolometer time-stream was clipped again at the 
3-cr level. 

All the datasets in the archive which were taken with 
the same chop throw were then made into a map using the 
SURF program REBIN. The maps were calibrated using jig- 
gle maps of both primary and secondary calibrators reduced 
in the same way. Figure 1 shows contour plots of the maps 
of these objects overlaid on images from the Digitised Sky 
Survey (DSS). The lowest level contour is la and increas- 
ing in la steps to 6a. Although it appears that the submm 
emission is concentrated in the centre of NGC5253, the dy- 
namical range of the optical data is much greater than the 
submm data, thereby making it difficult to draw any firm 
conclusions as to whether the optical emission is indeed more 
extended than the dust emission. 

The ffux was measured by using the DSS image as a 
guide to chose an aperture which contained as much of the 
galaxy and as little of the sky as possible. The flux errors 
were estimated using the procedure in Dunne et al. (2000). 
The fluxes and the apertures used are given in Table 1. 

Since the first draft of this paper another submm map 
of NGC1569 has been pubhshed (Lisenfeld et al. 2001). The 
new data was taken in much better conditions than the data 
we retrieved from the SCUBA archive and the flux measured 
is significantly higher than ours. We repeated our reduction 
of this object but the fiux determined remained 20% lower 
than the flux measured by Lisenfeld et al. (2001). The dis- 
crepancy may arise from the galaxy being larger than the 
SCUBA array, for which Lisenfeld et al (2001) have tried to 
correct. 



* The JCMT is operated by the Joint Astronomy Centre on be- 
half of the UK Particle Physics and Astronomy Research Coun- 
cil, the Netherlands Organization for scientific Research and the 
Canadian National Research Council 
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Table 1. SCUBA archive data and flux measurements 
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UM448 


72 


13 nov 1998 


4 


1280 


0.4 


34 


40 ±9 


4.1 


4.3 






13 nov 1998 


16 


1280 


0.4 










NGC5253 


3.2 


13 jan 1999 


77 


1280 


0.4 


41 


192±23 


30.5 


29.4 






13 jan 1999 


114 


1280 


0.4 










Mrk33 


24.9 


13 apr 1998 


15 


1280 


0.1 


36 


42±8 


4.7 


5.3 






13 apr 1998 


16 


1280 


0.1 










NGC4670 


11.0 


28 dec 1998 


62 


1920 


0.5 


38.6 


49±7 


2.6 


4.5 






21 jan 1999 


94 


1920 


0.4 














21 jan 1999 


95 


1920 


0.4 










NGC1569 


1.6 


10 sep 1997 


81 


1280 


0.4 


74 


269±32 


45.4 


47.3 






10 sep 1997 


82 


1280 


0.4 











Table 1. (1) Galaxy name (2) Distance to object from Nearby Galaxies Catalogue (TuUy 1988) (3) Date of observation (4) Telescope 
log run number (5) Integration time in seconds (6) tsso determined from skydips (7) Diameter of aperture used to measure 850/im flux 
(8) Measured 850/im flux (9) IRAS 60nm flux (10) IRAS lOO/im flux 



TABLE 2. Submm results 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Name 


D 


5850 




Nc 
N„ 


log Md(5860,K) 




(Mpc) 


(mJy) 


(K) 




{Mq) 


UM448 


72 


40 ±9 


45 


32.9 


7.41 


NGC5253 


3.2 


192±23 


45 


24.4 


5.39 


Mrk33 


24.9 


42±8 


42 


22.4 


6.51 


NGC1569 


1.6 


269±32 


40 


8.3 


4.91 


NGC4670 


11.0 


49±7 


38 


25.6 


5.87 


NGC2903 


6.3 


2120t 


34 


30.7 


7.02 


NGC4303 


15.2 


1180t 


41 


82.6 


7.54 


NGC6946 


5.5 


1200* 






7.70^ 


NGC7331 


14.3 


1900" 


34 


44.2 


7.69 


NGC1222 


32 


84 ± 16^ 


38 


6.0* 


7.02 


NGC7469 


64 


264 ± 30^ 


38 


12.3* 


8.13 


NGC7714 


37 


72 ± 13^ 


60 


13* 


7.05 


NGC3994 


41 


106 ± 20^ 


33 


9.8* 


7.34 


NGC3995 


43 


126 ± 24^ 


56 


498* 


7.47 


NGC5929 


33 


119 ± 30^ 






6.51 


NGC5953 


26 


182 ± 23'' 






7.20 


NGC5954 


26 


124 ±21^ 






7.03 


NGC6052 


62 


95 ± 15'' 


37 


16* 


7.65 


MCG+02-04-025 


122 


39 ± 18'' 


43 


6.2* 


7.83 


IR0335+15 


138 


44 ±9^ 


41 


14.5* 


8.01 


NGC5256 


109 


82 ± 17^ 


36 


8.6* 


8.07 


NGC7674 


113 


108 ± 20'' 


46 


87.7* 


8.25 



Table 2. (1) Galaxy name (2) Distance to object, determined using recession velocity taken from NED [The NASA/IPAC Extragalactic 
Database (NED) is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National 
Aeronautics and Space Administration] (3) 850/^m flux (f data from Amurc (priv. comm.), | data from Bianchi 2000, jj data from Bianchi 
1998, b data from Dunne ct al. (2000)) (4) Warm component dust temperature (5) Ratio of cold to warm dust, * results from Dunne 
& Ealcs (2001) (6) Dust mass calculated using equation 4, ^ dust mass taken from 200/im ISO data of Davies et al. (1999). NGC5929, 
5953 and 5954 are elements of interacting pairs that were not resolved by IRAS but were resolved by SCUBA. To calculate for these 
objects we fitted the Dunne & Eales two-component model to the 60, 100&850/im fluxes for both galaxies in the pair combined and then 
scaled the resultant dust mass using the 850/im fluxes for the individual galaxies. 



3 CALIBRATING e 

The fraction of metals incorporated in dust in the interstellar 
medium can be estimated from UV spectroscopic observa- 
tions of nearby stars, which show that a significant fraction 
of the metals that one would expect, from cosmic abun- 



dances, to be in the ISM are not there - they have been 
bound up in dust. If Ui is the number of atoms of a given 
chemical element relative to the number of hydrogen atoms, 
ai is the atomic weight of the element, and 5i is the propor- 
tion of that element bound up in dust, then 
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Figure 2. Two-component SED's assuming /3 = 2 for tlie IRAS 60 and lOO/xm and SCUBA 850/im fiuxes. Tlie solid line represents the 
composite two-component SED and the dot-dash lines show the warm and cold components. 



J2i ni X Si X ai 

J2i ^ «i 



(3) 



Recent UV observations of nearby stars with the Hubble 
Space Telescope have made it possible to calculate e more 
accurately because they have provided accurate depletion 
measurements for two of the most important elements, oxy- 



gen and carbon, strong evidence that the depletion is the 
same along different sight lines, and evidence that the cos- 
mic metal abundance is about two thirds the solar metal 
abundance (CardeUi et al. 1996; Meyer et al. 1998). To cal- 
culate e, we have assumed the solar metal abundance given 
in Pagel (1997) and then multiplied this by two thirds to 
obtain the cosmic metal abundances. It should be noted 
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TABLE 3. Dust and gas masses using metallicity dependent X-factors 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Name 


12 + log (O/H) 


Ref. 


log (Mhi) 


log (M^Ja 


log {Mh^)i 


log {Ma) A 


log (Md)/ 


















UM448 


8.08 


1 


9.40 


10.14 


10.92 


7.36 


8.08 


NGC5253 


8.34 


1,2 


8.29 


6.78 


7.07 


5.71 


5.73 


Mrk33 


8.40 


1 


8.47 


8.64 


8.94 


6.33 


6.54 


NGC1569 


8.16 


1,3 


7.75 


6.31 


6.92 


4.99 


5.04 


NGC4670 


8.30 


1 


8.22 


7.66 


8.10 


5.69 


5.83 


NGC2903 


9.12 


4 


8.98 


8.81 


8.03 


7.. 39 


7.21 


NGC4303 


9.01 


4 


9.42 


9.61 


8.99 


7.90 


7.64 


NGC6946 


9.06 


4 


9.49 


9.41 


8.72 


7.88 


7.69 


NGC7331 


9.03 


4 


9.66 


9.79 


9.14 


8.13 


7.87 


NGC1222 


8.57 


5 


9.08 






6.41 


6.72 


NGC7469 


8.80 


6 


9.18 


10.10 


9.80 


7.99 


7.76 


NGC7714 


8.53 


7 


8.93 


9.65 


9.75 


7.32 


7.41 


NGC3994 


8.61 


8 


9.45 






7.13 


7.13 


NGC3995 


8.66 


8 


9.79 






7.52 


7.52 


NGC5929 


8.18 


8 


8.63 


8.83 


9.46 


6.29 


6.77 


NGC5953 


8.73 


8 


8.76 


9.48 


9.29 


7.62 


7.20 


NGC5954 


9.16 


8 


8.63 


8.61 


7.77 


7.15 


7.18 


NGC6052 


8.65 


8 


9.58 


9.58 


9.71 


7.60 


7.67 


MCG+02-04-025 


8.87 


8 


9.24 






7.18 


7.18 


IR0335+15 


9.00 


8 




10.29 


9.69 


8.36 


7.76 


NGC5256 


8.75 


8 




10.46 


10.24 


8.28 


8.06 


NGC7674 


8.56 


8 


10.03 


10.67 


10.73 


8.39 


8.44 



Table 3. (1) Galaxy name (2) Oxygen abundance (3) References for oxygen abundances, 1. Lisenfeld et al. (1998) 2. Pagel et al. (1992) 
3. Skillman et al. (1989) 4. Zaritsky et al. (1994) 5. Petrosian et al. (1993) 6. Bonatto et al. (1990) 7. Gonzalez et al. (1995) 8. see § 4 (4) 
HI mass (5) H2 mass calculated using the X-factor of Arimoto et al. (1996) (6) H2 mass calculated using the X-factor of Israel (2000) 
(7) Dust mass estimated using eqn. 1 and the Arimoto et al. (1996) relation (8) Dust mass estimated using eqn. 1 and the Israel (2000) 

relation 



that the recent re-determination of the solar oxygen abun- 
dance by Allende Prieto et al. (2001), suggests that the solax 
metallicity is considerably closer to the ISM value than has 
been conventionally assumed. We have assumed the deple- 
tions given by Whittet (1992), except that we have preferred 
the recent measurements for carbon and oxygen obtained 
from the Hubble Space Telescope observations (Cardelli et 
al. 1996; Meyer et al. 1998). We obtain a value for e of 0.456. 



4 APPLYING THE METHOD 

Table 2 lists the submm results for all the galaxies that we 
are aware of which have metallicity measurements and sub- 
millimetre measurements. Many of the objects in the table 
are taken from the SCUBA Local Universe Galaxy Survey 
(SLUGS) of Dunne et al. (2000). Distances for objects with 
radial velocities v < 2000fcms^^ were taken from the Nearby 
Galaxies Catalogue (TuUy 1988). Table 3 gives the metal- 
licity measurements for the sample and also measurements 
of the masses of the atomic and molecular gas. Metallicities 
for some SLUGS objects were estimated using strong line ra- 
tios and the calibration method of Edmunds & Pagel (1984), 
relating R23 and [Nil]/ Ha and oxygen abundances in HIT 
regions. Line ratios were taken fron Keel et al. (1985) and 
Veilleux et al. (1995). Although more recent work has shown 
that the Edmunds & Pagel calibrator may over-estimate 
metallicities above solar, it is probably a reasonable esti- 
mator in the abundance range of these galaxies. 

There are two complications to applying the simple 



method we have outlined in the introduction: (a) the re- 
quirement of equation (2) for a temperature model; (b) the 
fact that estimates of gas mass depend to some extent on 
assumptions about metallicity. We will now discuss these in 
turn. 



4.1 Dust Masses 

Dunne & Eales (2001), in their investigation of the 
FIR/submm spectral energy distributions (SEDs) of galax- 
ies in the SCUBA Local Universe Galaxy Survey (SLUGS), 
found that the 450/850/xm flux ratio is remarkably constant 
for galaxies with a wide range of luminosities. They con- 
clude that the dust emissivity index (3 must be ~2. They 
then showed that with this emissivity index, there must be 
at least two dust components with different temperatures in 
galaxies, even in ultraluminous IRAS galaxies like Arp220. 

We fitted the Dunne & Eales two-component model 
to the five objects for which we have new SCUBA data. 
Since we only have three fluxes for each galaxy (60, 100 and 
850/im) , we fixed the temperature of the cold component to 
be 20K, the average of the values found by Dunne & Eales 
(2001). Figure 2 shows the fits for the five galaxies and Ta- 
ble 2 gives the temperature of the warm component and the 
ratio of the cold-to-hot dust (Nc/N^). We have also given 
the results of the two component fits for the SLUGS galaxies 
by Dunne & Eales (2001) and values for a few other galaxies 
for which we either made two-component fits or for which 
we took two-component fits from the literature. Prom the 
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(a) Dust mass on the ordinate calculated using the Ari- 
moto X-factor (eqn. 6) 




(b) Dust mass on the ordinate calculated using the Israel (c) Dust 

X-fa<;tor (eqn. 7) mass on the ordinate calculated using the Galactic 

'standard' X = 2.8 x 10^° H2cm-^/{kkms-^) 

Figure 3. Dust mass plots. The abscissa shows the dust mass calculated using equation 4 and the ordinate displays the dust mass 
determined with equation 1. The line is a least-squares best fit through the data. The stars represent the objects for which we reduced 
data from the SCUBA archive, the circles represent the SLUGS sample and the triangles represent galaxies whose fluxes were obtained 
from Amure (priv. comm.) & Bianchi (1998,2000). 



results of these fits we have estimated dust masses using 
equation (4): 



K850 



iVc 



Bs5oii',Tc) Bs5o{i',Tyj) 



(4) 



in which 5*850 is the flux density at 850/im, Bsso (i^, T) is 
the value of the Planck function at 850/im, Kgso is the 
mass-absorption coefficient at 850/im and D is the distance 
found for Qo = 1- We have initially adopted the value for 
Kd(850fim) of 0.077m^%'^ used by Dunne et al. (2000), al- 
though of course one of the goals of the present paper is to 
measure a new value for this. These initial estimates of the 
dust masses are given in Table 2. 



4.2 Gas Masses 

The mass of HI can be unambiguosly calculated from the 
intensity of the 21-cm line. Following Lisenfeld & Ferrara 
(1998), we have calculated the dust mass within the opti- 
cal disk, since the metallicity (and thus the dust content) 
outside the optical disk is likely to be low. 

Cold H2 is detected indirectly by means of a tracer such 
as CO. The CO — H2 conversion factor is however highly 
variable depending on temperature, density and metallicity. 
It is usual to determine the H2 mass in a galaxy by applying 
a 'standard' conversion factor, X = N{H2) / I(CO), where 
I(CO) is the velocity-integrated CO intensity and N{H2) is 
the molecular hydrogen column density. A complication here 
is that the conversion factor almost certainly depends on 
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-5.5 



-4 -3.5 -3 

(log (Md(850)/Mg) 



(a) Gas mass on abscissa calculated using the Arimoto X-fax;tor (eqn. 6). Plots using 
the Israel (2000) or Galactic standard X are almost identical. 



Figure 4. Dust-to-gas ratio and oxygen abundance. The stars represent the objects for which we took data from the SCUBA archive, 
the circles show the SLUGS (Dunne et al. 2000) sample and the triangles are the galaxies whose fluxes were obtained from Amure (priv. 
comm.) & Bianchi et al. (1998,2000). The filled squares are the data of Lisenfeld et al. (1998). The dot-dash lines indicate the distance 
that the Lisenfeld & Ferrara (1998) galaxies move on the plot when we use our submm fluxes rather than IRAS fluxes to estimate the 
dust masses. The two lines are the predictions of dust-formation models from Edmunds (2001), which have been normalised to the Galaxy 
(the open diamond). Both models imply that a constant fraction of metals are incorporated in dust at moderate and high metallicities. 
The solid line represents a model in which supernovae play an important role in dust production and the dashed line a model in which 
dust is only produced in some types of evolved low and intermediate mass stars. Some of the outlying points on these plots lack either 
HI or i?2 gas masses. With gas masses these data would move to tighten the correlation. 



metallicity. Wc investigated the dependence of our method 
on this uncertainty by trying a number of suggested rela- 
tionships between X and metallicity. 

Molecular hydrogen {H2) masses for the five objects 
from the SCUBA archive were calculated from the CO fluxes 
taken from the literature using 

(see Tinney et al. 1990) where Leo is in Kkms^^p(? . Ico 
values were taken from Sage, Salzer, Loose and Henkel 
(1992) and from Taylor, Kobulnicky and Skillman (1998). 
Leo was calculated from the relation Leo = {'Tr^)Ico 
where r is the radius of the telescope beam on the galaxy 
in pc. Gas masses for the SLUGS sample were taken from 
Dunne et al. (2000). For those objects whose 850/im fluxes 
were obtained from Amure (priv. comm.) and Bianchi et al. 
(1998, 2000) we took Sco values from Young et al. (1995). 

We investigated the dependence of the method on the 
metallicity-dcpcndcncc of the X-factor by trying two dif- 
ferent suggested relationships for the dependence of X on 
metallicity. Arimoto, Sofue and Tsujimoto (1996) propose 
that X has a strong dependence on metallicity. 



logX' = -1.0(12 -I- log[0]/[H]) + 9.30 (6) 

where X' = NhJIco x H2 / {Kkms-^) . 

Israel (2000) notes the probable dependence on ultra- 
violet ratiation field, but also suggests a simplified approxi- 
mate metallicity dependence for X, 

logX = 12.2- 2Mog[0]/[H] (7) 

We have used both of these relations to determine H2 masses 
for our sample. These are listed in Table 3. As an ad- 
ditional test of the sensitivity of the method to the X- 
factor we have also used a "standard" conversion factor, 
X = 2.8 X H2cm-'^ l[kkms~^], from Bloemen et al. 
(1986) derived from observations of the Galaxy. Table 3 also 
gives the dust masses calculated using equation 1. 



5 INTERPRETATION &: DISCUSSION 

In figure 3 we have plotted the dust masses calculated by the 
two different methods against each other. We first performed 
a simple least-squares fit on the data points. All of the least 
squares solutions were consistent with a slope of unity. This 
suggests that the properties of dust, in particular Ka and e, 
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are similar for dwarf galaxies and large galaxies like our own. 
On the assumption that this is true, the value of Kd can be 
determined in the following way. If the value of Kd we have 
initially assumed for plotting Fig. 3 is correct, the points 
should lie on a lino of slope unity passing through the origin. 
By measuring the actual offset of the points from the line, 
we can calculate a new value of Kd- We obtain values of 0.07 
± 0.02, 0.07±0.01 and 0.07±0.01 rn^kg'^ using Arimoto, 
Israel and Galactic CO — H2 X-factors respectively. 

The method clearly does not depend critically on the 
metallicity-dependence of the CO-to-i/2 X-factor, presum- 
ably because either the overall contribution of H2 to the 
gas-mass is in all cases relatively small, or, where the molec- 
ular content is high, the X-factors do not vary greatly at 
these metallicities. The main uncertainty in the method is 
the basic assumption that the fraction of metals in dust (e) 
is a universal constant. If this is not correct, our estimate of 
this constant, from observations of the interstellar medium 
in our own galaxy will not be applicable to other galaxies. 
However, Figure 3 is itself evidence that this is not too bad 
an assumption. The scatter around the best fit lines is only 
a factor of 2 in the ordinate. This scatter includes all ob- 
servational errors, as well as galaxy-to-galaxy variations in 
Kd and £. On the extreme assumption that all the scatter 
is caused by variation in e, we estimate that the galaxy- 
to-galaxy variation in e can only be a factor of ~ 2. If this 
extreme assumption is correct, then our estimate of k, would 
have an error of a factor of ~ 2 for an individual galaxy. 

Our conclusion that the fraction of metals incorporated 
in dust is the same for giant and dwarf galaxies is in contra- 
diction with the conclusions of Lisenfeld & Fcrrara (1998), 
who found that metals are less effectively incorporated in 
dust in dwarfs than in giant galaxies. Figure 4 shows metal- 
licity plotted against dust-to-gas mass ratio for our sample 
and that of Lisenfeld & Ferrara. The Lisenfeld & Ferrara 
points suggest a non-linear relationship between dust-to-gas 
ratio and metallicity. However we have re-observed some of 
the objects from their sample with SCUBA and find that 
their dust-to-gas mass ratios are increased when estimated 
from submm rather than far-infrared flux. This is expected 
since in their study, Lisenfeld ct al. (1998) used the IRAS 
60/100/im flux ratio to determine a dust temperature and 
hence may have underestimated the dust mass by missing 
cold dust at longer wavelengths. The rc-obscrved points axe 
consistent with a slope of unity. Very recently Lisenfeld et 
al. (2001) have suggested that a low dust mass is appropri- 
ate for NGC1569, whether such a model is applicable to all 
dwarf galaxies, or is unique to NGC1569 is not yet clear. We 
have also plotted two of the galaxy dust evolution models of 
Edmunds (2001) in figure 4. The solid line shows a model in 
which supernovae play an important role in dust production 
and the dashed line a model in which dust is produced in 
some kinds of evolved low and intermediate mass stars. Our 
current dataset is unable to make a distinction between these 
two models but further submm observations at low metallici- 
ties {12 + log{0/H) < 7.8) should help demonstrate whether 
or not condensation of dust grain cores in supernova ejecta 
is an important dust source. Such observations will be quite 
challenging as the number of objects suitable for observation 
at 850/im which have very low metallicities is small. 

In Figure 5 we have reproduced a plot from Alton (2001) 
that shows estimates of the dust mass absorption coefii- 
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Figure 5. The dust mass absorption coefRcient in the submm (re- 
produced from Alton ct al. (2001)). The data is from Hildebrand 
1983, Casey 1991, Rengarajan 1984, Sopka ct al. 1985, Alton et 
al. 2000, Agladzc et al. 1994, Draine & Lee 1984, Boulanger et al. 
1996, Bianchi et al. 1999 and this work. The data from Agladze 
et al. 1994 is for amorphous dust grains at 20K. The line shows 
K <x u'^ 

cient in the IR-Submm from a variety of sources. Sopka et 
al. (1985) made observations of the thermal emission from 
the dust envelopes of evolved stars. Rengarajan (1984) re- 
ports results from observations of centrally heated infrared 
sources deeply imbedded in molecular clouds. Boulanger et 
al. (1996) studied the dust-to-gas correlation at high Galac- 
tic latitude using COBE data. The Draine & Lee (1984) 
points come from a model for diffuse Galactic dust. Both 
the Hildebrand (1983) and Casey (1991) points are from 
Galactic reflection nebulae. The Alton (2000) measurement 
is from SCUBA observations of the spiral galaxy NGC891 
while Bianchi ct al. (1999) used COBE & IRAS data. Fi- 
nally, Agladze et al. (1994) performed laboratory experi- 
ments on forsterite {Mg2SiOi). The line shows k. oc z/^. Our 
data (open triangle) is consistent with the COBE data of 
Boulanger et al. (1996) as well as that of Draine & Lee 
(1984), Hildebrand (1983) and Bianchi et al. (1999). 



6 CONCLUSIONS 

We have presented a new method for determining the 850/Lim 
dust mass absorption coefficient, and obtain a value of 
0.07±0.02 m^kg-'^. 

Our method for determining this coefficient avoids the 
problems encountered by previous investigators (Hildebrand 
1983, Bianchi 1999, Alton 2000). Using SCUBA observations 
at 850/im we can probe the dust emission directly. We do not 
make any eissumptions about the properties of the dust nor 
do we rely on similarity between dust in Galactic reflection 
nebulae and that in the general ISM. Our only assumption 
(admittedly an important one) is that the fraction of metals 
incorporated in dust is a universal constant. This seems a 
justified assumption however since the small scatter in Fig. 
3 suggests that the fraction of metals in dust is the same for 
a wide range, from dwarf to high mass galaxies. 

Having used both metallicity-dependent (Israel 2000; 
Arimoto et al. 1996) and independent (Bloemen et al. 1986) 
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CO — H2 conversion factors to derive a molecular gas mass 
we have concluded that the metaJlicity-dependence of the 
CO — H2 X-factor is not an important consideration in this 
derivation of Kd- 

The COBE data of Boulanger ct al. (1996) are con- 
sistent with our value for ksso as are the data of Draine & 
Lee, Bianchi and Hildebrand. The /tsBo value of Alton (2000) 
that was derived from SCUBA observations of the edge on 
galaxy NGC891 is ~ 3 times higher than ours. They do how- 
ever admit that uncertainty in their technique means that 
their quoted value for the dust emissivity could be in error 
by a factor of 2 either way. Agladaze et al. (1994) laboratory 
experiments yield a Kgso that is ~ 2 times higher than ours, 
however they use (3 = 1.3 for wavelength extrapolations. 
Sopka et al. 1985 also use /3 = 1.2. 

Unfortunately, the metallicity range of our sample does 
not extend low enough to make a distinction between the 
dust formation models proposed by Edmunds (2001). We 
intend to acquire 850/xm data for the rest of the Lisenfeld & 
Ferrara (1998) sample where possible. 
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